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Abstract 

We derive semiclassical ground state solutions that correspond to the quantum 
Hall states earlier found in the Maxwell-Chern-Simons matrix theory. They realize 
the Jain composite-fermion construction and their density is piecewise constant as that 
of phenomenological wave functions. These results support the matrix theory as a 
possible effective theory of the fractional Hall effect. A crucial role is played by the 
constraint limiting the degeneracy of matrix states: we find its explicit gauge invariant 
form and clarify its physical interpretation. 



1 Introduction 



In this paper we continue the study of noncommutative and matrix gauge theories [1][2][3] 
that can be used as effective non-relativistic theories* [4] [5] of the fractional Hall effect [6]. 
The original Susskind proposal was based on the noncommutative Chern-Simons theory 
[7]: it developed in the Polychronakos Chern-Simons matrix model [8], that was analysed 
by several authors [9] [10] [11] [12] [13] [14] [15][16]. This model is a U(N) gauge theory 
in one dimension (i.e. time) with Chern-Simons kinetic term; it involves two Hermitean 
matrix coordinates, X\(t),X2(t), that are noncommuting, [Xi,^] = i B I, where is a 
constant "background charge" and / is the identity matrix. In the semiclassical limit [7], 
the matrix model describes two-dimensional incompressible fluids in strong magnetic field 
B with Laughlin's values of the filling fraction, v = l/(k+ 1), where k = B0 is integer 
quantized [17]. 

At the full quantum level, some problems were found in matching the matrix model 
to the physics of the fractional Hall effect: one point was that the theory reduced to the 
eigenvalues, X a , a = 1, . . . ,N, of X = X\ + i Xi, does not really describe electrons in the 
lowest Landau level with coordinates X a ; owing to matrix noncommutativity, the Laughlin 
wave function is deformed at short distances [11][15]. Another issue was that this model 
cannot easily describe the general Hall states with Jain's filling fractions: V = mj (mk± 1), 
m = 2,3,... [18]. 

In our earlier paper [19], we showed that these problems can be overcome by upgrading 
the Chern-Simons model to the Maxwell-Chern-Simons matrix theory. This includes an 
additional kinetic term quadratic in time derivatives and the potential, V = — gTr[Xi,X2] , 
with the coupling g > that controls matrix noncommutativity. All terms in the action 
are determined by the gauge principle because they can be obtained by dimensional reduc- 
tion of Maxwell-Chern-Simons gauge theory, as discussed in the literature of DO-branes in 
string theory [3]. 

The Maxwell-Chern-Simons matrix theory reduces to the earlier Chern-Simons model 
for large values of B (with finite g). However, in the different g = 00 limit (with finite 
B), corresponding to [Xi,^] = 0, it does provide a sensible physical description of the 
fractional Hall effect: after reduction to eigenvalues, one finds electrons in Landau levels 
interacting with a two-dimensional potential of Calogero type, l/\X a — Xb\ 2 (see section 
2.2). For general B < °o, the additional coupling g in the theory allows one to interpolate 
between matrix (g = 0) and electron (g = <*>) dynamics. 

Furthermore, the theory is exactly solvable at g = [19]: it was found to describe a 
*We refer to [15] for an introduction to noncommutative theories of the quantum Hall effect. 
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matrix extension of the Landau levels, where the gauge-invariant many-body states are 
given by matrix generalizations of Slater determinants. Although the degeneracy of matrix 
states grows exponentially with energy, it was possible to control it by introducing suitable 
projections in the theory. We showed that the constraint (A a t,) m ps 0, \fa,b, projects the 
theory to the "lowest m matrix Landau levels", with m = 1,2,. . .; the m-th reduced theory 
naturally possesses non-degenerate homogeneous ground states with filling fractions of 
the Jain series, v = m/(mk+ 1) [19] (see section 2). Indeed, the solutions of the gauge 
invariance conditions and of the constraint A m ps give raise to a rather surprising structure 
of ground states that corresponds to the Jain composite-fermion construction of ansatz wave 
functions [18]. 

These g = states exactly match the phenomenological Jain wave functions under matrix 
diagonalization, that is formally achieved at g = °°. Therefore, we conjectured that these 
matrix ground states have smooth deformations for g > into the physical g = 00 states, 
namely that no phase transitions are found for < g < °° when the system is at the specific 
densities [19]. 

A number of problems remained to be further investigated: 

• Understand the matrix Jain states, e.g. compute their densities and observables quan- 
tities. 

• Understand the projections (A a t,) m ps in more physical terms. 

• Study the phase diagram of the theory for < g < °° at the relevant densities. 

In the present paper, after reminding earlier results [19] (section 2), we find the gauge 
invariant form of the projection A m ps and its semiclassical physical meaning in terms of 
single-particle occupancy (section 3). Next, we study the matrix Jain states in the semi- 
classical approximation, by analytically solving the classical equations of motion, further 
constrained by the Gauss law and the semiclassical version of the A m ps condition (section 
4). The ground states are found to be two-step droplets of incompressible fluid with piece- 
wise constant density; this is the same density shape of the phenomenological Jain states 
before projection to the lowest Landau level [18] (where the density of incompressible 
fluids becomes strictly constant). 

The fact that the matrix Jain states at g = already have the expected droplet density of 
physical g = °° states, supports our earlier claim that these ground states could remain stable 
while varying < g < °° [19]. Other ground states corresponding to generalized Jain con- 
structions with different filling fraction, although possible in the g = theory, are found not 
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to possess piecewise constant density. We argue that the modulated density shape is a sig- 
nal of ground- state instability at finite g values, since the corresponding phenomenological 
Jain states (g = °°) are known to be unstable [18]. We complete our study of semiclassical 
solutions by describing the quasi-holes excitations above the matrix Jain states. Finally, 
in the conclusion we briefly discuss the ways to study the Maxwell-Chern-Simons matrix 
theory for g > 0. 



2 Jain states in Maxwell-Chern-Simons matrix gauge the- 
ory 

2.1 Lagrangian and Hamiltonian 

In this section we recall the matrix theory of quantum Hall states proposed in [19] and the 
derivation of ground states that are matrix analogs of the Jain composite-fermion states. 
The theory involves three time-dependent N x N Hermitean matrices, Xi(t), i = 1,2 and 
Ao(t), and an auxiliary complex vector The Lagrangian contains a Maxwell-Chern- 
Simons kinetic term, a uniform "charge background" and the \\f "boundary term" of 
ref.[8]: 



/ 



dtTr 



in "i B 2 o 

-{D t Xif + -EijXiDtXj + 8 -[X u X 2 } 2 + B0A o 



-iJ^D t y. (2.1) 

The covariant derivatives are: D t X t = X t - i [A ,X ; ] and D t \\f = ?A y. Under the U (N) 
gauge transformations: X, — > UXiU', Aq — > U (Aq — id/dt) W, and \\r — > U\\r, the ac- 
tion changes by a total derivative, such that invariance under large gauge transformations 
requires the quantization, B0 = k E 7h [17]. Hereafter we set m = 1 and measure dimen- 
sionful constants accordingly. 

The variation of S w.r.t. the non-dynamical field Ao gives the Gauss-law constraint; its 
expression in term of coordinates X t and conjugate momenta n ; , i = 1,2, reads: 

G w , G = i [Xi,ni] + i [X 2 ,U 2 ] - B0 I + y®y f . (2.2) 

The trace of G fixes the norm of the auxiliary vector \|/, 

TrG = — > ||y|| 2 = B0/V = &V. (2.3) 



We note that \\r has trivial dynamics, \\r(t) = V|/(0) = const., and it is necessary to represent 
the Gauss law on finite-dimensional matrices that have traceless commutators [8]. In a 
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gauge in which all \\f components vanish but the last one, the term (kl — \|/®\|/^) in (2.2) 
becomes the "traceless identity", k L7, L7 = diag( 1 , • • • , 1 , 1 — N) . 

In gauge-invariant quantization, all 2N 2 matrix degrees of freedom X' ab axe quantized 
and the Gauss law is imposed on states: G generates U (N) gauge transformations of X, 
and \|/, and G = implies that physical states are U (N) singlets subjected to the additional 
condition (2.3) fixing the total number of \|/ a components equal to kN. The Hamiltonian 
can be written: 

H = B Tr (a 1 " a) + *N 2 - |Tr [X u X 2 ] 2 , (2.4) 
after introducing the variable, 



A = ^(x 1+ ix 2 ) + ^(n 1+ m 2 ) 



(2.5) 



and its adjoint A''', involving the magnetic length, i = a/2/B. These quantities obey the 
commutation relations of A^ 2 independent harmonic oscillators: 



A a b-,A 



cd 



ad Vbc 



[[AabAcd]] = 



(2.6) 



Therefore, for g = the Hamiltonian describes Landau levels populated by A^ 2 two-dimensional 
"particles" with phase-space coordinates, {JV^X^}, a,b = 1, . . . ,N, i= 1,2. Degenerate 
angular momentum excitations are described by an independent set of oscillators: 



B = ±(x 1 -ix 2 ) + ^(n!-/n 2 ) , 



abi D cd 



= §ad§bc , [[Bab, B cd\] = 0; 

that commute with A, A" 1 ". The total angular momentum of the A^ 2 particles is 
J = Tr(Xi n 2 - X 2 ni) = Tr (b^b - A f A) . 



(2.7) 
(2.8) 

(2.9) 



The N 2 - particle states are further constrained by the Gauss law (2.2): as we shall see 
later, this enforces a kind of generalized exclusion statistics. It is convenient to introduce 
the complex matrices, 



X = X\ + i X 2 , X = X\ — i X 2 , 

n = I(n!-m 2 ), n = ^(ni + in 2 ), 



(2.10) 



and use the bar for denoting the Hermitean conjugate of classical matrices, keeping the 
dagger for the quantum adjoint. Hereafter we set the magnetic length to one, i.e. B = 2. 
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Figure 1 : Phase diagram of the Maxwell-Chern-Simons matrix theory. The vertical axes 
g = and g = oo have been discussed in ref. [19]. The Chern-Simons matrix model [8] is 
found at B — > °o in the left down corner. 

For states with constant density''', the angular momentum measures the extension of the 
"droplet of fluid", such that we can relate it to the semiclassical filling fraction v c i by the 
formula, 

In a physical system of finite size, one can control the density of the droplet by adding a 
confining potential Vc to the Hamiltonian: 

H -> H + V c = H + ©Tr^fl) + co„ Tr (fi +n B n ^j . (2.12) 

The strength (0 is of order 0(B/N) such that the structure of Landau levels is not destroyed 
by putting he = 0(N) particles per level. The higher order terms 0((D n ) also commute 
with the g = Hamiltonian: one can show that their eigenvalues on constant-density states 
are very large 0(N n ~ l ) for fillings hb > n and thus they can be used to simulate finite-box 
boundary conditions [19]. 

In Figure (1) we illustrate the phase diagram of the Maxwell-Chern-Simons matrix the- 
ory as a function of the parameters B/m and g. The analysis of [19] found the properties 
of the theory on the axes g = and g = °°. For g = 0, the theory is exactly solvable and 
possess a set of ground states that are in one-to-one relation with the Laughlin and Jain 
ground states with filling fractions v = mj (mk+ 1). These states are selected by the value 
of k and by adjusting the density via the parameters (0, (0„ of the confining potential (2.12); 



^See section 4 for the definition of density in matrix theories. 
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furthermore, they are unique non-degenerate states obeying the conditions (A a b) mx H = 
that projects out degeneracies specific of matrix states, as it will be explained momentarily. 

The limit B — > °° at g = (and g finite) in the Hamiltonian (2.4) corresponds to vanishing 
kinetic term, i.e. to the constraints A = A* = of lowest Landau level [20]. In this limit, 
the Gauss law (2.2) reduces to: 

G = -iB [Xi,X 2 ]-Be + Y<g>V' , (2.13) 

and it uniquely fixes the noncommutativity of matrices. The potential term in H (2.4) 
becomes a constant for all physical states and the theory reduces to the Chern-Simons 
matrix model [8] with trivial dynamics, H = Vc [19]. 

For g = °o (B finite), the theory describes ./V ordinary electrons with coordinates given 
by the complex eigenvalues A, a , a = 1, . . . ,N, of X, interacting with the two-dimensional 
Calogero potential. 



2.2 g — > oo limit and electron theory 

For large g values, the potential term in H (2.4), g Tr[X,X] 2 , restricts the configuration 
space of complex matrices to those commuting with their conjugate, the so-called "nor- 
mal" matrices [21]. Therefore, X can be made diagonal by unitary (gauge) transformation 
and the theory completely reduces to the eigenvalues, that are interpreted as electron coor- 
dinates. Any complex matrix can be written [22] as, X = U (A + R)U, in terms of a unitary 
matrix U (the gauge degrees of freedom), a diagonal matrix A (the eigenvalues) and a upper 
triangular complex matrix R (additional d.o.f.). Therefore, the gauge-invariant degrees of 
freedom different from the eigenvalues contained in the R matrix are suppressed for g — > oo. 

We can thus take the diagonal gauge for X and decompose the momenta II, II in diagonal 
and off-diagonal matrices, respectively called p and T: 

x = a, n = p+r, TT = p+r. (2.14) 

In this gauge, the Gauss law (2.2) can be solved at the classical level and it determines the 
off-diagonal momenta, 

Upon inserting them into the Hamiltonian (2.4), the diagonal and off-diagonal components 
decouple and one obtains, 

■t-#-*)M^iCT' <2 - l6) 
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The quantization can be done on the remaining independent variables, which are the com- 
plex eigenvalues X a and their conjugate momenta p a [23]. Therefore, the theory reduced to 
the eigenvalues corresponds to the ordinary Landau problem of N electrons plus an induced 
two-dimensional Calogero interaction. The measure of integration on matrices (2.18) also 
reduces to that of ordinary electrons after incorporating one Vandermonde factor in the 
wave functions [21]. This causes a renormalization of the filling fraction from the semi- 
classical expression (2.11): l/v = 1 + 1/Vci 

We conclude that the Maxwell-Chern-Simons matrix theory in the g = 00 limit makes 
contact with the physical problem of the fractional quantum Hall effect: the only difference 
is that the Coulomb repulsion e 2 / r is replaced by the Calogero interaction k 2 / r 2 . Numerical 
results [24] [25] [18] [26] indicate that quantum Hall incompressible fluid states are rather 
robust and do not depend on the detailed form of the repulsive potential at short distance, 
for large B, at least for the qualitative features. The g = 00 theory is not, by itself, less 
difficult than the ab-initio quantum Hall problem: the gap is non-perturbative and there are 
no small parameters. The advantage of embedding the problem into the matrix theory is 
that of making contact with the solvable g = limit. 

2.3 Matrix Jain states at g = 

The wave functions of the Maxwell-Chern-Simons theory take the form, 

Vj/ = e -Tr(r*)/2-w/2 <&(X,X,Y) , (2.17) 

and their integration measure reads: 

(¥i|¥ 2 ) = J 'DX'DX 2)\|/©\j7e _ ™-W^)*(x,x,\|/) <£ 2 (X,X,V) • (2.18) 

At g = 0, the energy and momentum eigenstates are obtained by applying powers of the 
A\ b and B\ b operators (2.5,2.7) to the empty ground state = exp(-TrXX/2-\jA|//2) 
(as in ordinary Landau levels), leading to: 

vj/ = e -Tr xx /2- w/2 <$>(B,A,\\f) , E = BN A , j =N B -N A . (2.1 9) 

The wavefunction 4> is a homogeneous polynomial of B = X — d/dX and A = X — d/dX, 
that can be treated as onumber matrices because they commute among themselves. The 
energy E = BNa and momentum J = Nb — Na of the state are expressed in terms of the 
eigenvalues of Tr A^A and Tr B^B, Na and Nb respectively, that count the total number of 
matrices A and B in 4>. 
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In the lowest Landau level, the wave function obeys A a ^ = 0, V a,b, thus the polynomial 
part 4> does not contain any A ab : it is a analytic function of the B ab variables, equal to the 
X a t,. The physical states 4>(X, \|/) obeying the Gauss law are U (N) singlets that contain any 
number of X a t and Nk copies of the \\f vector (owing to (2.3)). The solutions for k = 1 are 
given by [10]: 

*{„!,...,*,} = Z ai - aN (X n ^) ai ---(X n ^) aN , 0< ni <n 2 <---<n N , (2.20) 

where £ is the completely antisymmetric tensor and {«,•} any ordered set of integers. So- 
lutions for k > 1 are obtained by multiplying k terms (2.20), leading to the expressions, 
„i \ /„* Ji \- The ground state in the confining potential TriXX^) is given by the 

\n l ,...,n lf f"\n 1 ,...,n N } <-' ^ ^ 

closest packing {0, 1,. . . ,iV — 1} that has the lowest angular momentum, i.e. lowest degree 
inX: 



(2.21) 



If we diagonalize the complex matrix X by similarity transformation: X = V~ l A V, 
A = diag{X\, . . . , Aw), \|/ = V -1 (|), the dependence on V and (|) factorizes and the powers of 
eigenvalues make up the Vandermonde determinant A(k) = Y[ a <b(^a — ^b), as follows: 

<D fe , g ,(A,V, V ) = (detV)- fc n (K~h) k (U <$>c) ■ (2.22) 

\<a<b<N V c J 

The central piece is indeed the Laughlin wave function [24], upon interpreting the eigen- 
values as electron coordinates [10] [11]. The value of the filling fraction (2.11) is: 

v = FTT ' (2 ' 23) 

by keeping into account the extra Vandermonde coming from the integration measure. 

Eq.(2.22) is the most interesting result obtained in the noncommutative approach and 
the Chern-Simons matrix model [8] [10]: that of deriving the Laughlin wave function from 
gauge invariance in a matrix theory with background charge 9. Furthermore, Susskind's 
semiclassical analysis [7] showed that, in the limit 0^0, this matrix state describes an 
incompressible fluid in high magnetic fields with density, p = 1/(2710) , and classical 
filling fraction v c i = 2np /B = 1 /B0 = l/k in agreement with the earlier identification. 

The analog of the Laughlin quasi-holes are realized by shifting the occupation numbers 
of matrices, e.g. by the state in (2.20) with {ni,n2, ■ ■ ■ ,«m} = {1>2, • • • ,N}. This has Aj? = 
0(N) and thus a finite gap AE = 0(B). On the other hand, the quasi-particle excitations 
cannot be realized in the Chern-Simons matrix model [8]. In general, states with higher 
density do not exist in this theory, because they would need to populate higher Landau 
levels that are absent. 
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(a) (b) (c) 

Figure 2: Graphical representation of gauge invariant states: (a) general states in the lowest 
Landau level (cf. Eq.(2.20)); (b) and (c), N = 3 examples in the second and third levels 
involving both matrices, B (thin line) and A (in bold). 

Further difficulties in matching the Chern-Simons matrix model to the Laughlin physics 
at the quantized level were discussed in Refs. [11][15][13]: since the matrices X,X are 
noncommuting (cf. (2.13)), the theory is not well suited for the description of electron 
degrees of freedom. The reduction [15] to eigenvalues (X fl ,X fl ), interpreted as electron co- 
ordinates, shows that the measure of integration (2.18) does not become that of the Landau 
levels and the Laughlin state (2.22) gets deformed at short distance. These findings should 
be contrasted with the results in section 2.2 for Maxwell-Chern-Simons matrix theory, that 
possesses a well-definite physical limit at g = °o (although difficult to solve). 

The states (2.20) can be represented graphically as "bushes", as shown in Fig. (2a). The 
matrices B a t (i.e. X a j,) are depicted as oriented segments with indices at their ends and 
index summation amounts to joining segments into lines, as customary in gauge theories. 
The lines are the "stems" of the bush ending with one \\r a , represented by an open dot, and 
the epsilon tensor is the N-vertex located at the root of the bush. Bushes have N stems of 
different lengths: n\ < ri2 < ■ ■ ■ < n^. The position if of one B on the £-th stem, 1 < it < tig, 
is called the "height" on the stem. 

The solutions of the Gauss law (2.2) for states in higher Landau levels are $ polynomials 
involving both B and A matrices: given that they transform in the same way under the gauge 
group, these polynomials can be represented by bushes whose stems are arbitrary words of 
B and A, as shown in Fig. (2b, 2c), where A a t, is represented by a bold segment. Since 
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two stems cannot be equal, one obtains a kind of Fermi sea of N "one-particle states" 
corresponding to the N strands. However, there are additional degeneracies with respect to 
an ordinary fermionic system, because in each stem all possible words of A and B of given 
length yield independent states (for large N), owing to matrix noncommutativity, as seen in 
fig.2(b) and 2(c). 

The complete filling of all the available degenerate E > states clearly gives very dense 
and inhomogeneous fluids that are incompatible with the physics of the quantum Hall ef- 
fect. Moreover, the matrix degeneracies lead to a density of states in the g = theory that 
grows exponentially with the energy. This is a characteristic of string theories that is not 
suitable for the Hall effect [3]. Of course, for g > the potential Tr[X,X] 2 restricts the 
matrix noncommutativity and reduces the degeneracy: at g = °°, this is not present and the 
theory can describe a physical electron system, as is clear from the discussion of section 
2.2. 

Given that the g > theory is very difficult to solve, in ref. [19] we introduced a set 
of projections that limit the matrix degeneracy at g = and are explicitly solvable. These 
projections are expressed by the following constraints on the wave function, 

-=-J 4>(A,fl,y) =0, Va,fc, (2.24) 

for a given value of m. The m = 1 case is the lowest Landau level discussed before with 
no A dependence, while m taking successive values m = 2,3,... gradually allow larger 
A multiplicities and thus matrix degeneracies. Note that in equation (2.24), each matrix 
component A ah is raised to the m-th power, without index summation: the condition is 
nevertheless gauge invariant and it admits an equivalent manifestly invariant form that is 
discussed in section 3. 

The results of ref.[19] were rather interesting: not only the projections (2.24) allow 
homogeneous ground states suitable for describing quantum Hall fluids, but also they pre- 
cisely occur in the Jain pattern of filling fractions, v = m/(mk+ 1), and their derivation 
repeats step-by step the Jain "composite fermion" construction [18]. 

Let us recall the main points of the analysis of [19]. Consider first the projection (2.24) 
for m = 2: the solutions are polynomials that are at most linear in each component A a j,. 
Let us imagine that one or more A matrices are present at points on the bush as in Fig. (2). 
The differential operator (2.24) acts by sequentially erasing pairs of bold lines on the bush, 
each time detaching two strands and leaving four free extrema with indices fixed to either 
a or b, with no summation on them. For example, when acting on a pair of A located on 
the same stem, it yields a non-vanishing result: this limits the bushes to have one A per 
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stem at most. The A 2 pa conditions can be satisfied if cancellations occur for pairs of A 
on different stems, as it follows: 

(a*) <&=••• + ^••'•••'•••(•••M /(( /V /( ,---V // MV' / ') + ••• , (affixed). (2.25) 

This expression vanishes forM = N due to the antisymmetry of the epsilon tensor. The 
general solution of (2.24) is given by bushes involving one A per stem at most (max N 
matrices in total), with all of them located at the same height on the stems [19] . In formulas: 

^ {nu ... m+u ..., nM} = e'i-'" n ft (B P AB"^) ik , 

k=i k=e+i 
< n\ < • • • < ri£ , < ri£ + i < • ■ ■ < ■ (2.26) 

If the matrices A, B were diagonal, these states would be Slater determinants of ordi- 
nary Landau levels. The matrix states have further degeneracies by commuting A, B pairs: 
however, commutations are almost impossible in the the solution (2.26) of the A 2 pa con- 
straint, were it not for the p dependence. This shows how the projection works in reducing 
matrix degeneracies. 

The ground state in the A 2 pa theory with finite-box conditions corresponds to homo- 
geneous filling all the allowed states in the first and second Landau levels with N/2 "gauge 
invariant particles" each. It reads: 

N/2 N/2 

*i/2 > * = e' 1 " iw ft [B k V. II P^"V). , (2-27) 

k=\ lk k=\ V 7 l N/1+k 

with angular momentum J7 =N(N — 4)/4. This state is non-degenerate due to the vanishing 
of the p parameter in (2.26). It has filling fraction v* = 2, assuming homogeneity of its 
density*. 

The ground states for k > 1 are products of k bushes: they obey the constraint A 2 ps 
provided that the two derivatives always vanish when distributed over the bushes. Given 
one bush of type (2.27), obeying A 2 gs = 0> one can form the state, 

^Jfc+l/2, gs = ®k-l, gs &1/2, gs , (2-28) 

where the other (k—l) bushes satisfy A ^-1, gs = and actually are Laughlin's one (2.20). 
The angular momentum value for this state corresponds to the filling fraction 1/v = k+ 1/2 
(cf(2.11)). 

We thus find the important result that the A 2 pa projected Maxwell-Chern-Simons the- 
ory possesses non-degenerate ground states that are the matrix analogues of the Jain states 



*This will be shown in section 4. 
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obtained by composite-fermion transformation at v* = 2, i.e. 1/v = 1 /v* + k. The matrix 
states (2.28,2.27) would actually be exactly equal to Jain's wave functions, if the A, B ma- 
trices were diagonal: the \\f dependence would factorize and the matrix states reduce to the 
Slater determinants of Jain's wave functions (before their projection to the lowest Landau 
level) [18] [26]. 

The correspondence extends to the whole Jain series: the other v* = m non-degenerate 
ground states, corresponding to V = m/ (mk+ 1), are respectively obtained in the theories 
with A m m projections. They read: 

®k+l/m, gs = ®k-l, gs ®l/m, gs > (2.29) 

where, 



N/m 

q _ ch-iN FT 

l/m, gs c 11 

k=\ 



' h ^ ' h+N/m ^ 7 l k+(m-l)N/m_ 



(2.30) 

Note that in the A m ps theory, the lower density states that were non-degenerate in the 
A^ ps theories, k < m, become degenerate; actually, when the boundary potential is tuned 
for letting (2.30) to be the ground state, the former states (k < m) become excited states. In 
conclusion, in ref.[19] we found that the ground states of the properly projected Maxwell- 
Chern-Simons matrix theory reproduce the Jain pattern of the composite fermion construc- 
tion [18]; the matrix states are non-degenerate for specific values of the density that are 
controlled by the boundary potential. 

These results indicate that the Jain composite fermions have some relations with the 
DO-brane degrees of freedom and their underlying gauge invariance. Both of them have 
been described as dipoles: according to Jain [18] and Haldane-Pasquier [27], the compos- 
ite fermion can be considered as the bound state of an electron and a hole (a vortex in the 
electron fluid). On the other side, matrix gauge theories, and the equivalent noncommuta- 
tive theories [3], describe DO branes that are point- like dipoles in the low-energy limit of 
string theory. Although further relations between these two dipole pictures remain to be 
found, in this paper we shall find further evidences that the matrix ground states describe 
semiclassical incompressible fluids with some of the properties of Hall states. 

2.4 Generalized Jain's hierarchical states 

In the A m projected theories with m > 3, there are other solutions of the Gauss law for 
k > 1 besides the Jain states (2.29). These are obtained by combining products of any k = 1 
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m 


#/>/ = 1 


#Pi = 2 #pi = 3 #pi = 4 


l/v = l + l/v c , 


E/B 


1 


k 


k+1 





2 


k-l 


1 


k+1/2 


N/2 


3 


k-2 


2 


k 


N 


3 


k-l 


1 


k+1/3 


N 


4 


k-3 


3 


Jk— 1/2 


3N/2 


4 


k-2 


1 1 


k-1/6 


3N/2 


4 


k-l 


1 


£+1/4 


3N/2 



Table 1: Examples of standard (2.29) and generalized (2.31) Jain states for fixed value 
of k, ordered by the level m of projection, A m ps 0, and their filling fraction v and energy 
E (disregarding the confining potential). Note that the experimentally relevant values are 
* = 2,4. 

solution (2.30), as follows: 




A^i gs = 0, 9=1 + E(pi-1). (2-31) 

(=l 

In this equation, we also wrote the associated filling fractions using Eq.(2.1 1), i.e. assuming 
homogeneous densities, and the condition A q that they satisfy. The Jain mapping to 
a single set of V* = q effective Landau levels does not hold for these generalized states 
(having more than one pi > 1). Actually, analogous generalized states were considered by 
Jain in his composite-fermion theory [18], but were discarded due to their instability (small 
or vanishing numerical gap). 

Let us compare the generalized (2.31) and standard (2.29) Jain states at fixed values of 
the background k (keeping in mind that the physical values are k = 2, 4). The energy of the 
generalized states is additive in the V* = pu k = 1, blocks and reads: 

£_L + ... + _L,,, = ^ £(W-1)+V C . (2.32) 

The analysis of some examples shows that these additional solutions have in general higher 
energies for the same filling or are more compact for the same energy than the standard Jain 
states (2.29) (see Table 1). States of higher energies are clearly irrelevant at low temper- 
atures; higher-density states strongly deviate from the semiclassical incompressible fluid 
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value v = l/(k+ 1) for background B9 = k, that is specific of the Laughlin factors [7]. The 
analysis of the corresponding semiclassical states in section 4 will show that most of these 
states are not incompressible fluids (density not piecewise constant). 



3 Properties of the projection A m « 

In this section we discuss the physical meaning of the projection: 

(A ab ) m *¥(A,B) = 0, Va,b, (3.1) 

that limits the degeneracy of matrix quantum states at g = 0. Although the operator (A a b) m 
is not gauge invariant, its kernel restricted to gauge invariant states yields a gauge invariant 
condition^, as explicitly seen in section 2. Therefore, there should exist a manifestly gauge 
invariant expression for this condition, that is found in this section. 

A simple example is useful to clarify the following discussion. In a two dimensional 
quantum mechanical problem with rotation invariance (0(2) global symmetry), we con- 
sider the condition: 

(-) \ m 
-\ 4>(r 2 ) =0, r 2 =x 2 +y\ (3.2) 

where $ is a reduced (polynomial) wave function. The condition is not 0(2) invariant 
but its kernel acting on rotation invariant functions does: indeed, it limits the order of the 
polynomial to C^r" 1-1 ). This example suggests two remarks: 

• The condition (3.2) can have many different forms, that correspond to points on its 
orbit in the "gauge" 0(2) group: for example, an equivalent form is (d/dy) m <& = 0, 
corresponding to a 7T./2 rotation. All these conditions are equally satisfied. 

• A manifestly gauge-invariant expression can be obtained by integrating over the 
gauge orbit, as follows: 

P m — > P% 1 - = J Q dd (cos0— + sin0 ^;) • (3-3) 

However, this vanish for m odd: the average looses information because the operator 
(d/dx) m is not positive definite. Clearly, it can be made positive (and gauge invariant) by 
contracting with another gauge-dependent term to obtain powers of the dilatation operator 

D'" = (x i d/dx i ) m . 

8 A formal proof of this statement is given in Appendix A 
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We are now going to follow analogous steps for the condition A m 0. First we find 
an equivalent, more general form. Consider an infinitesimal SU (N) gauge transformation 
U = 1 + ieT: the Hermitean matrix T can be expressed by the matrices 2s with a single 
non- vanishing component, = § l a 5 J b , in symmetric or antisymmetric combinations, T = 
2s +2?( Ji ) or T = i(E^ — E^). Upon performing the gauge transformation, the m = 2 
constraint (2.24), (U*AU) 2 ab , acquires an additional 0(e) term that should also vanish on 
the wave functions obeying, (A ab ) 2 m = 0: 

» A ab [E'J,A] ab = A ab (8 ai A jb -A ai 8 jb ) , Vi^j, V a, b. (3.4) 

We now analyse the various cases: 

• I. If a = b and i = a or j = a, we obtain the conditions, 

~ A aa A ja w A aa A ai , ViJ^a. 

• II. If a 7^ b, we obtain, 

1 . for i = a and j 7^ b — > ps A ab Aj b , V 7 ^ a, , 

2. for / 7^ a and 7 = & — > A a £ A aj - , V i^fl,i, 

3. for / = a and j = b — > p=j A ai (A to - A aa ) . 

Note that each term in the linear combination of case II. 3 independently vanishes by 
case I. 

These conditions can be summarized as follows: 

AabK'b^ = 0, Va,a',b, 

A ab A aU ^ = 0, Va,b,b'. (3.5) 

They are more general than the original expression (2.24) for m = 2, corresponding to 
a = a' or b = b'. Of course, iteration to 0(e 2 ) of the gauge transformation produce further 
identities: these involve linear combinations of A 2 terms that are not particularly useful; for 
example, one such condition is: A ab Aj c + A^ac ps 0. 

The 0(e 2 ) analysis is necessary to obtain the generalized constraint for m = 3: the 0(e) 
expression is similarly, A ab A ab A a y ps 0, and its further transformation yields, 

« 2A ab [E^^^Aay +A ab A ah [E ij ,A] a y ■ 
This expression contains the m = 3 constraint analogous to (3.5): 

A a b Aa'b A a"b ^ = 0, Vc,fl',/,fe, 

A ab A M A aV ,^ = 0, Va,b,b',b", (3.6) 
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together with other relations involving linear combinations of cubic terms. 

Following the 0(2) example, we can now transform the new expressions (3.5) into posi- 
tive definite operators. We recall that the lowest Landau level condition corresponds to the 
vanishing of the total energy, that is a sum of positive terms: 

H = Tr(A* a) = Y, A *ab A ab ~ & A ab ^0,Va,b. (3.7) 

a,b 

We can construct the following positive definite expressions: 

& = £ Al a A\ a A ab A aV , (3.8) 

a,b,b' 

Ql = E A L A la A ab A a'b , (3.9) 
a,a',b 

whose vanishing is equivalent to the m = 2 conditions (3.5). These quantities are not yet 
gauge invariant but are convenient for the physical interpretation. We introduce the (gauge 
variant) energy operators for one-particle matrix states, that are summed over matrix indices 
of one row or column of A ab , Z a or Z' b , respectively: 

Z a = £ A la A ab , 4 = £ A l A ab • (3-10) 

b a 

Using these energy operators, we can rewrite (3.8,3.9) as follows: 

q 2 = Y^z a (z a - i) , ^ = E4(4-i)- (3- 11 ) 

a b 

In this form, the constraints Qfl! = = admit the following physical interpretation: 
there is a gauge choice in which the allowed states contains at most one "particle" in the 
second Landau level (energy equal to one) for (a,b) indices belonging to each row and 
column. 

The constraint for m = 3 (2.24) similarly becomes: 
Q 3 = £ Z a (Z a - 1) {Z a - 2) , <2' 3 = £ Z' b (Z' b - 1) (4 - 2) . (3.12) 

a b 

This requires that there at most 2 particles in the second Landau level or a single particle 
in the third level for any set of indices in a row or column. The matrix labels are not gauge 
invariant, then these occupancies are only verified in specific gauges; nevertheless, the 
present form of the constraints can be implemented in the semiclassical limit on expectation 
values, {Aab), as explained in section 4. 

Next, we obtain the gauge-invariant form of the constraint Q2, Q' 2 by averaging over the 
gauge group. We define: 

Gf = (<DU Q' 2 {U) =£/*©[/ U\. Al uljA] a A ak U kh A dl U lb , (3.13) 

b J 
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where T)U is the invariant Haar measure [28]. The integrand is positive definite for any 
U value, because it can be thought of as the norm of a vector: QziU) ~ Y*b \A • v^| 4 , 
where vi^ = U an $n * s a rotated unit vector. Therefore, we do not loose any information by 
performing the group average. 

Group integrals of products of U,U^ matrices can be found e.g. in ref.[28]: their re- 
sults can be described as follows. Representing the unitary matrices with upper and lower 
indices, U a t, — > U b , [U^) a b — > [U^) a b , the result of integrating n (U,U^) pairs is a combi- 
nation of n delta functions relating the upper indices among themselves times other n deltas 
connecting the lower indices. The simplest integral is: 

f <du (u'y a , u b , b = U ab d a , h> . 

In the general case of n {U , U^) pairs, the pairings of upper (lower) indices by delta func- 
tions follow patterns given by the permutation of n elements, with specific weight for each 
conjugacy class of permutations [28]. For n = 2, one finds: 

d ab d cd d a>b ,d c , d , + d ad d c %, d ,d c , b , - 1 (5 fl *5"W5 c ^ + 5^5^,5^) . 

In the case of the constraint Q' 2 (3.9), all the upper indices are simultaneously taking the 
same value b; thus, the different delta-function pairings of upper indices take the same unit 
value. As a result, the pairings of lower indices get averaged over, and reduce to a plain 
sum over all pair permutations: 

C ~ (hi hj + hj hi) A] a , A) a A ak A a ,{ . (3.14) 

Upon commuting the operators to bring summed indices close each other, we finally find 
the manifestly gauge-invariant form of the A 2 ks constraint (disregarding the normaliza- 
tion): 

2f'~0, Q% L = Tr (aUaU) + (jrA^y - (N+ 1) Tr (a + a) . (3.15) 

The same expression is also obtained by group averaging the other operator Q2 in (3.8). 
One can check that the action of the gauge-invariant constraint Q 8 2 u on bush wave functions 
(cf. section 2) is completely equivalent to that of the gauge- variant condition A 2 « [19]. 

The gauge invariant form of the m = 3 constraint can be similarly obtained by group 
averaging (3.6), leading to: 

<t = I Kb A \y A i»> K (2) y \ m b ■ o.iq 



1 
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The form of this expression corresponding to (3.15) is not particularly illuminating. The 
gauge-invariant expression (3.16) straightforwardly generalizes to higher m values. 

In conclusion, in this section we have found equivalent forms of the projections A m « of 
g = matrix states: the first expression (3. 1 1 ,3- 12) in terms of occupation numbers is useful 
for the semiclassical limit considered in the next section; the second expression (3. 15,3. 16) 
is manifestly gauge invariant. In the latter form, the constraint can be added to the Hamil- 
tonian with a large positive coupling constant to realize a softer form of projection, where 
matrix states violating the constraint are now allowed but possess very high energy. For 
example, the quasi-particles excitations over the Jain ground states V = m/ (mk+ 1) would 
be possible, while they are absent in the A m O-projected theory, as explained in section 
3. A detailed analysis of this issue is postponed to a following publication. 

4 Droplet ground state solutions 

In this section we study the g = Maxwell-Chern-Simons theory in the semiclassical limit: 
we solve the classical equation of motion including the quantum constraints, first for the 
ground states and then for the quasi-hole excited states. We shall find the semiclassical 
states that correspond to the quantum states with homogeneous filling and composite- 
fermion structure of section 2 [19]. The motivations for this semiclassical analysis are 
twofold: on one side, previous experience [8] [3] [12] [13] [16] with noncommutative field 
theory has shown that the classical fluid dynamics incorporates some properties of the full 
quantum theory. From another side, it is know that the Laughlin states in the quantum 
Hall effect are incompressible fluids that become semiclassical in the thermodynamic limit 
N — > oo [29] [30]. The semiclassical ground states we find in this section are also incom- 
pressible fluids which, we believe, may give rather accurate descriptions of the quantum 
matrix states for large N values [21]. 

Let us start by writing the classical equations of motion: the Hamiltonian of the Maxwell- 
Chern-Simons theory at g = can be written as follows: 

H = 2Tr(AA) + roTr(flfl) + Tr [A ([A,A] + [£,£]- + 

+ £r fl (Z fl -y) + ^T' b (z' a -i) , Y ,Y = 0,l,...,m- 1(4.1) 

a b 

We set B = 2, B0 = k E Z, and included the Gauss law constraint via the Hermitean La- 
grange multiplier A. The projection A m ps analyzed in section 3 is enforced by adding two 
other Lagrange multipliers r fl , T' b times the energies, Z a = £ fc A ba A ab , Z' b = E« Ha Kb, 
of single-particle states with matrix indices summed over rows or columns. We replace 



18 



the nonlinear constraints (3.11,3-12) with linear expressions involving the parameters y,Y 
taking the allowed values of Z a ,Z' b . Since the constraints are not gauge invariant, we shall 
assume that we work in a gauge where they take integer values. The gauge-invariant form 
of the constraint (3.15) found at the end of sections 3 is also not convenient because it 
would lead to non-linear equations of motion that cannot be solved analytically. For the 
same reason, we limit the confining potential (2.12) to the quadratic term: later we shall 
see how to avoid ground state degeneracies that may arise with this potential. 

We vary the Hamiltonian with respect to A, B, canonically equivalent to the original X , n, 
and obtain the equations: 

ik ab = 2A ab - [A,A] ab + A ab (T a + ry , (4.2) 

iB = -[A,B] +oo5 , (4.3) 

G = [A,A] + [B,B] -k + \\f(g>W = , (4.4) 

Z a = £ A ba A ab = y, y=0,l,...,ra-l . (4.5) 

b 

4 = Y, I baA ab =y f , Y = 0,l,...,m-1 . (4.6) 

a 

We first discuss ground state solutions corresponding to A = B = 0. 



4.1 Laughlin ground states 



We can gain some intuition on the matrix equations (4.2,4.3,4.4,4.5) by recalling the solu- 
tions of the Chern-Simons matrix model [8], corresponding to the lowest Landau level in 
our theory. In the case A = A = 0, we have B = X and the ground state equations reduce to: 



[A,B] = oo 5 , 



G= [B,B] -& + \|/®\j7 = 



(4.7) 



These are the commutation relations of a truncated quantum harmonic oscillator, with A 
playing the role of Hamiltonian. We can diagonalize it by gauge choice and write standard 
oscillator matrices in the "energy" basis, |0), |1 ),..., |iV— 1) [8]: 

N-l N-l 

B = Vk ^y/n\n)(n — 1|, A = 03 ^n\n)(n\ , 

n=l 

X|/ = v^v | iV— 1> . 



n=l 



(4.8) 



In matrix form for /V = 3: 





( 





0\ 




/ 1 





0\ 




(0\ 


B = \fk 


1 








, A = 03 





2 





, \|/= VkN 









V2 






\o 





3 ) 




V 1 / 



(4.9) 
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The solution is characterized by the angular momentum, J = Tr B B = k N(N — l)/2 and 
by vanishing energy. 

A good definition of the density of semiclassical fluids in matrix models [8] is given in 
terms of the gauge invariant eigenvalues of R 2 = XX, 

N-l 

P(r 2 ) = £ 5(r 2 -o ; ) , c ( GSpec(tf 2 ) . (4.10) 

In the limit N — > °°, this becomes a piecewise continuous function that describes two- 
dimensional rotation-invariant distributions (p(r) = p(r 2 )/7t). A discrete approximation 
is: 

P(^ = E^V,„ (4.1D 

involving the Kronecker delta and the ordered set of distinct a, eigenvalues, O; < a,, i < j, 
with multiplicities 

In Polychronakos' solution of the Chern-Simons matrix model (4.8), the matrix R 2 is 
already diagonal, 

R 2 =XX =diag(0,*,2*,...,(#-l)Jfc) . (4.12) 

Its density (Fig. 3) is constant and describes a circular droplet of fluid with Laughlin filling 
fraction [7]: 

V - = ^ = F wi,h p " = 2^e' (4 - 13) 

Note that an ordering ambiguity in the definition of R 2 was resolved by matching to the 
angular momentum spectrum. As said before, there is a shift in the filling fraction in 
quantum theory due to a contribution from the integration measure, 1/v = 1 + l/v c /; thus, 
k should be even for describing electrons. 



4.2 Jain ground states 

As we showed in section 2.3, the Maxwell Chern-Simons theory contains Jain-like ground 
states (2.29), that involve higher Landau levels (A 0). Their filling fractions can be written 
as in composite fermion construction [18], 

1 1 

- = — + k+l, k even, v*=2,3,..., (4.14) 

and their energy and angular momentum values are recalled in Table 1. We first note 
that these states are characterized by energies 0(N) and angular momentum J = 0(N 2 ), 
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Figure 3 : Plot of the density for the semiclassical Laughlin ground state (4.8) for v c / = l/k, 
k = 4, and N = 400. 

thus implying that the matrix A must have elements of 0(1) and be much smaller than B. 
Indeed, the constraints, Z a ,Z' b = 0, 1, . . . ,ra — 1, limits the squares of A ab matrix elements 
summed over each row or column to take at most the total value (m — 1). Were it not for this 
constraint, the A, B matrices could be rescaled in the ground state equations (4.2,4.3,4.4) to 
eliminate the k dependence, leading to solutions with E = 0(kN) at least. 

We now describe the solution of the ground state equations of motion in the A 2 
projected theory (y,Y = 0, 1 in (4.5)). Under some minor hypotheses, we find a single 
solution corresponding to the unique quantum state with v* = 2 (2.28) [19]. Working in 
analogy with the Laughlin case (4.9), we shall try a distribution of R 2 eigenvalues leading 
to a piecewise constant density. We can consider the gauge in which A is diagonal, A = 
diag(£ a ), and assume that V|/ has a single non-vanishing component, i.e. the last one, as in 
(4.9), such that the term, (k I — \|/<g)\j7), is also diagonal. The equation for B (4.3), 

(£ a -i b )B ab = (oB ab , (4.15) 

requires that B is a raising operator, i.e. non vanishing on a single diagonal, B ab oc 8 a ^ + „; 
moreover, the Gauss law (4.4) requires, [B,B] ~ k I, apart from 0(1) corrections due to 
[A, A]. Therefore, B should be non- vanishing on the first diagonal: 

B ac = 5 a , c+ i b c+ i , c = 0,...,N-2. (4.16) 

Eq. (4.15) implies evenly spaced A eigenvalues, £ a +i — £ a = G>, and leaves the components 
b c undetermined. The equation (4.2) for A reads: 

A ab ± — T a + T' b = (a-b)(0 - 2 , (4.17) 

that can always be solved for T a ,T' b . The constraints, Z a ,Z' b = 0, 1, imply that A a t has 
one non-vanishing element per row and column, at most, equal to one. If it had exactly 
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one element per row and column, it would be the representation of a permutation, a E Sn- 
Therefore, we can write: 

A ab = §a,c(b) a b+\ ■> <* b = 0,l , a e S N ■ (4.18) 

We now consider the Gauss law (4.4): all terms in this equation are diagonal matrices, 
leading to a system of (N — 1) scalar equations for the A, B matrix elements {a.b,bb\. Note 
that both matrices AA and BB are diagonal and thus their elements are positive integers in 
the semiclassical theory: b\ e Z+. After introducing, 

fo = b 2 b , a b = a 2 b , (4.19) 

we obtain the system: 

Pi = fc-OCi+Otyi), 

P2-P1 = fc-a 2 + a o(2 ), 

• • • ) 

= £-ocjv-i +oc o(A ,_i) • (4.20) 

The solution can be found by thinking to the expected shape of the droplet. The quantum 
state (2.28) is made of k generalized Slater determinants with homogeneous filling of ./V 
"one-particle" states'". Each one-particle state is expected to give a constant contribution 
to the density of the droplet: there are (k— 1) Laughlin terms and one term with N/2 
"particles" in the second Landau levels spanning half of the angular momentum range, as 
confirmed by the quantum numbers, E = BN/2 and J = (k - 1 + 1 /2) N 2 /2 + 0(N) . The 
contribution are additive in terms of angular momentum eigenvalues, J ~ Tr BB = Y,f=i P( 
(the 0(N) contribution of Tr AA is subdominant for Af — > °°). Therefore, we expect, (3 { - ~ 
(k— l)i , for one half of the range, say < i < N/2, and (3/ ~ (k+ l)i for the other half. 
Moreover, p\ should be continuous at i = N/2 in order to obey the corresponding equation 
witha, = 0(1). We take: 

p. = (k+ 1) i — N, ^<i<N. (4.21) 

The solution for A is found by inspection: it has N/2 non- vanishing elements equal to one 
on the diagonal of the lower half sector. 

""This Fock-space analogy is meaningful for diagonal matrices, and may not be correct in general: its 
limitations will be discussed in section 4.3. 
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Summarizing, the ansatz semiclassical ground state solution for v* = 2 is given by (N 
even): 



5 



A = 



N/2 N-l 

£ y/n(k-l) | n){n- \ \ + £ v / n(k+ 1) -N | n) (n - 1 



A' 



E 

n=0 



n + -)(n 



(4.22) 



In matrix form for N = 4, it reads: 



5 



( 











o\ 




/o 








o\ 


































v/2(*-l) 










1 











\ 








V3&-1 


0/ 




\o 


1 





0/ 



(4.23) 



This solution has same energy £ = BN/2 of the quantum state (2.28) and same angu- 
lar momentum J = (k— \/2)N 2 /2 + O(N) to leading order (cf. Table 1). The matrix 
R 2 = (B +A) (B +A) contains off-diagonal terms from the mixed products: however, these 
give subdominant 0(1/ VN) corrections to the eigenvalues as is clear in a simple two-by- 
two matrix example. Thus, Spec(7? 2 ) = Spec (55) (1 +0(l/y/N)), confirming the earlier 
identification of droplet shape with angular momentum spectrum. 

In Fig. 4(a), the density (4.11) of the droplet of fluid is plotted by computing the exact 
spectrum for Af = 400: up to finite-iV fluctuations, this is a two-step constant density as 
anticipated. We recall that the same droplet shape is found for the Jain phenomenological 
states before their projection to the lowest Landau level [18]; the density becomes constant 
only after projection 1 1 . 

In ref.[19], we argued that the matrix ground states at g = match one-to-one the phe- 
nomenological Jain states that are good ansatz in the physical limit g = °°: the two sets of 
states become identical in the limit of both X,X diagonal, that can be formally reached at 
g = oo. To establish a relation at the quantum level, we would need to consider the evo- 
lution of the matrix ground states as the coupling is varied in between, < g < °°, and to 
check that the gap never vanishes, i.e. that there are no phase transitions in (5,g) plane 
(cf. Fig.l) separating the g = and g = 00 regions at these density (i.e. total angular mo- 
mentum) values [19]. While this behaviour remains to be proved, it is supported by the 
result that the matrix (g = 0) and phenomenological (g = oo) states have similar densities 
of incompressible fluids. 



II The lowest-level projection in the matrix theory cannot be done at present, lacking an understanding of 
the g > regime: at g = 0, it would give a trivial result because the Laughlin state is the unique lowest-level 
ground state for any k value. 
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Figure 4: Plot of the density for the Jain matrix ground states with l/v c / = 1/v* +k, for 
k = 4 and N = 400: (a) v* = 2 (4.22); (b) v* = 3 (4.24); and (c) v* = 4 (4.26). 
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We also note that the solution (4.22) could also be obtained in the lowest-level theory 
(Chern-Simons matrix model) by replacing the A matrix with N /2 different "boundary" 
auxiliary fields \\f — > \|/ a , a= 1, . . . ,N/2. This multi-boundary generalization of Polychron- 
akos' model has been considered in ref.[9]: it naturally describes multicomponent droplets, 
i.e. 1/v = n/kfov n boundary fields. However, the description of Jain states is rather unnat- 
ural, because the number of auxiliary fields is macroscopic and should be adjusted for each 
Jain state; moreover, this theory does not admit the physical limit of commuting matrices. 

The solution (4.22) can be easily generalized for the theory with projection A 3 0, 
possessing a Jain ground state with v* = 3: this is found at the specific density that is 
reached by tuning the boundary potential. The constraint now allows the A a t, components 
(4.18) to take values a h = 0, 1, \fl\ we assume again a single non- vanishing element per 
row and column, eq. (4.18), otherwise the commutator, [A, A], would have off-diagonal 
terms that cannot be matched in the Gauss law equation (4.4). Therefore, the equations 
(4.20) are unchanged. Let us recall that the quantum solution contains (k— 1) Laughlin 
terms and the v* = 3 piece that puts three particles in the same angular momentum state, 
ranging from zero to N/3. Thus, the B ansatz contains eigenvalues spaced by (k— 1) for 
2/3 of the droplet and by (k+2) for 1/3 of it. The matrixA that solves the Gauss law (4.20) 
involves elements on a diagonal extending for 2/3 of the matrix (N should be a multiple of 
3). In conclusion: 



2N/3 N-l 

B = £ y/n(k-l) \ n){n- \ \ + £ y/n(k + 2) -2N\n)(n-l \ , 

IN 
3 



n=\ 



£ \n+-)(n\ + £ V2\n + -)(n 

n=0 J N 5 

3 



(4.24) 



In matrix form for Af = 6: 
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(4.25) 



The droplet shape plotted in Fig.4(b) has again two steps, up to local fluctuations that vanish 
for Af — > oo. 
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The ansatz solution with v* = 4 in the theory A 4 xs again involves a matrix B with two- 
speed spectrum and a matrix A with elements a b = 1 , V2, \/3, on the diagonal extending 
over 3/4 of the matrix (iV multiple of 4): 



3N/4 



N-l 



B = £ y/n(k-l) | n)(n-l | + £ v / "(^ + 3 ) -3# !")("- 1 I , 



w=l 



--1 



2JV 



-1 



3,V 



-1 



A = £ + £ v^|n + ^)<n| + £ v^|n + j)(n| . (4.26) 



n=0 

In matrix form for N = 8: 
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(4.21) 



The density for N = 400 and k = 4 is plotted in Fig.4(c). 



4.3 Correspondence of semiclassical and quantum states 

Here we provide a simple argument to support the identification of the semiclassical so- 
lutions with the quantum states of section 2. Consider first the correspondence for the 
Laughlin states, (2.21) and (4.8). We choose the gauge in which the expectation values of 
B matrix elements on the quantum state take the classical values (4.8) found in the previous 
section, up to subleading corrections for N — > °o. Let us rewrite the N = 4 wave function in 
terms of these non- vanishing terms only**: 

®k, gs = 

~ [e 3210 y 3 (5 23 V3) (5i 2 5 2 3X)/3) (B iB 12 B23^3)] k ■ (4.28) 

** Although this expansion should hold for N — > °°, we write the N = 4 case for simplicity; the expression 
for general N can be easily inferred. 
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This "semiclassical wave function" describes "particles" with matrix indices, (01), (12), (23), 
in angular momentum states that precisely match the occupation numbers B ab B ba given by 
the classical solution (4.12), equal to (k,2k,3k), respectively. This is a self-consistent ar- 
gument for the correspondence of states: in the semiclassical N — > °° limit, the quantum 
states match the semiclassical solutions for the leading occupation numbers. 

A similar relation holds for the v* = 2,3,4 Jain states. For v* = 2, the quantum wave 
function is (2.28); we evaluate it on the semiclassical non- vanishing A a b,B a b values (4.22), 
given explicitly for N = 4: 



®k+l/2, gs 



=2 s r=3 



k-l 



e a ia2 a 3 a 4 ^ (g>y)^ 
xe a ia2 a 3 a 4 ^ (J^ (JB^)^ 

~ [£ 3210 ^3 (5 23 V3) (5i 2 5 23 V3) {BoiBn^^)] 

X£ 3210 y 3 (5 2 3^3) (Ai 3 \|/ 3 ) (Ao 2 5 2 3V3) • (4.29) 

The "one-particle" occupancies of both energy and angular momentum states given by the 
wave function again match the expectation values of the corresponding number operators, 
AabAba and B a t,Bb a , of the classical solution. The correspondence extends to the other v* = 
m states that have spectrum of occupancies given by (4.24,4.26). This argument support 
our belief that the large ./V limit of the matrix theory is semiclassical for the incompressible 
fluid ground states (piecewise constant density) and their small excitations. 



4.4 Generalized Jain states 

In the analysis of [19], we found other quantum solutions to the constraint A m m 0, for 
m > 3, besides Jain composite fermion wave functions. They were recalled in section 2.4, 
eq. (2.31) and summarized in Table 1: these are analogs of Jain's generalized hierarchical 
states, made by products of two or more wave functions with higher-level fillings (pi > 1). 
In the semiclassical analysis, we find that some of these states have corresponding solutions 
with piecewise constant density, while most of them do not. Besides, we find spurious 
ground states that are allowed by the simplistic quadratic boundary potential used in (4.1). 
Let us describe these solutions in turn. 



4.4.1 Spurious solutions 

There is a variant of the composite-fermion solution for m = 3,4, Eqs. (4.22,4.24), 
where the A matrix elements take the same values, but their positions are permuted. For 
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m = 3, this is: 

N/3 N-l 
B = £ ^n{k-2) | n)(n-l \ + £ y / n(k+ \) -N \ n) (n- I | , 

A = £ v^|„+^)(„| + £ \ n + ^)(n\ . (4.30) 

The total energy and angular momentum values are the same as those of the m = 3 Jain state, 
E = BN, 3 ~ (k — 2/3)N 2 /2 (cf. Table 1). The corresponding B matrix again describes 
a two-step droplet. In order to find the corresponding quantum state, we use the classical- 
quantum correspondence of the previous section. The single-particle occupation numbers 
of the classical solution, for e.g. N = 6, correspond to those of (k — 2) Laughlin fac- 
tors plus the occupations (3,6,9), (2,2, 1, 1) and (12) for the components (523,^34,545), 
(A02A13A24A35) and \|/ 5 , respectively. These components should fit into two wave func- 
tion of the type (2.26) that obey A 2 (the product wave function obeying A 3 0). The 
solution, rewritten in gauge invariant form, is: 

4> = (<*>!, ^ (>1«2«3«4«5«6 ^ (g 3^ (Atf 2 ^ (^Vk) 

x (e™*^ ^ ( 5 ¥)fl2 ( A ¥ ) fl3 (AB» a4 (A5 2 X|/)a 5 (A5 3 ¥)a 6 ) • (4-3D 

In this state, we recognize that some strands do not have minimal length: thus, this is an 
excited state for a Hamiltonian with more realistic boundary terms (2.12) realizing finite- 
box conditions''"''. The expectation value of the higher boundary potential (Tr (b 2 B 2 j ) on 
this state is actually larger than that of the m = 3 Jain state (4.24) with same energy and 
angular momentum: this confirms our interpretation of the solution (4.30). 

4.4.2 Other two-step density states 

Among the generalized Jain state in Table 1, there are those made by two kinds of terms, 
as follows: 

% gs = (*!,„)*- (*!,„)", I = £ + (*-») + !, n = 2,3,..., 

(4.32) 

in the notation of Eq. (2. 31). They obey, A n+l rs 0, forn = 2, 3, . . ., and violate the composite- 
fermion transformation (4.14) [18]. In the droplet interpretation of classical solutions of 
section 4.2, we expect a density of R 2 eigenvalues equal to (k + n) for half of the spectrum 
and (k — n) for the second half. The ansatz has the two-block structure of the solution 
(4.22), with maximal value A a t, = n in agreement with the constraint. 

^The quadratic potential used in (4.1) is known to yield such degeneracies ([29]). 
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Figure 5: Plot of the density for the generalized Jain states: (a) 1 /v c / = l/k—1 (4.33); (b) 
l/ Vcl = l/k-3/2 (4.35), with k = 4 and N = 400. 



The first non-trivial value is n = 2, i.e. m = 3 in Table 1: 



N/2 



N-l 



B = £ y/(k-2)n\n){n-l \ + £ y/(k + 2)n-2N \n)(n- 1 



n=\ 



-1 



£ V2|^ + n)(/i 

rc=0 



In matrix representation for N = 4: 
B = 



o 



^2(^-2) 



A = 



y / 3k^2 







t/2 





The analogous state for n = 3, corresponding to (^m ^y)" 3 (^l, 3 , is 



N/2 



N-l 



B = £ y/(k-3)n\n)(n-l \ + £ y/(k + 3)n-3N \ n)(n- 1 



«=l 



£ v^|^ + n)<n|, 



and in matrix form for N = 4: 



B 




A = 







V3 

V3 



(4.33) 



(4.34) 



(4.35) 



(4.36) 



o o v^/t-i) o 

In Fig. 5, we plot the density for these generalized Jain states: these are droplets with two- 
step constant density similar to that of composite-fermion states. At present, we do not have 
strong arguments to dispose of these additional states: this issue will be further discussed 
in the conclusion. 
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r 2 

Figure 6: Plot of the density of the 4- step excited state in the A 4 ps theory, \/v c l ~ 
k- 1-0.14, with k = 4 and N = 400. 

4.4.3 States with many-step density 

Other generalized matrix Jain states (2.31) in Table 1, for A m ps 0, m > 4, are made 
by the product of three or more different terms. The simplest one for m = 4 is 4> v = 
4>*~ 2 <J>!/ 2 4>j/ 3 with energy E/B = 3N/2 and angular momentum 7 ~ (& - 1 - 1 /6)N 2 /2. 
Within the droplet interpretation of classical solutions discussed before, we seek for a three- 
step solution (N multiple of 6), 

N N 2N IN 

P,-~i(*-2), 1< / < — ; p,-~i*, y< z '<^-; p f ~i(*+3), —<i<N. 

However, there is no simple A a b solution with entries (0, 1, y/2, \/3), that fulfills the Gauss 
law equation for the same energy and angular momentum of the quantum state. We proved 
this fact for an ansatz with piecewise constant density making up to 6 steps. A four-step 
solution (see Fig. 6) can be found with quantum numbers differing macroscopically from 
the quantum values, E ~ 1.4 N,J ~ (k — 1 — 0.14)iV 2 /2: in particular, the larger angular 
momentum identifies it as an excited state. Presumably, the quantum state can be better 
approximated by allowing a very large number of steps, leading to a modulated (or singular) 
density profile in the large-./V limit. This result indicates that most of the multi-component 
generalized matrix quantum states found in [19] for projections A m ks 0, m > 4, are not 
semiclassical incompressible fluids. 

4.5 Quasi-holes solutions 

As recalled in section 2.3, the g = matrix theory, projected by A m ps 0, possess quasi-hole 
excitations above the v* = m Jain ground states. In this section we give the corresponding 
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semiclassical solutions for v* = 2, corresponding to deformation of the density of solution 
(4.22) in fig.4(a). 

The classical equation of motion for A and B, eq. (4.2,4.3), are linear and admit a general 
solution for excitations: 

Aabif) = e-^ r " +2 ^ (e ltA A(0)c- ltA ) e~< r »> , 

V / ab 

B(t) = e itA 5(0) &- itA . (4.37) 

Therefore, we should only solve the Gauss law (4.4) and the constraint (4.5). 

In the two-step fluid density in fig.4(a), one can have more than one quasi-hole corre- 
sponding to punching either of the two possible fluids. A hole in the complete fluid is 
obtained by generalizing the quasi-hole of the Laughlin state found in ref.[8]: it is a de- 
formation of the B matrix (4.22) that describes a quasi-hole of charge q > situated at the 
origin. The matrix A remains unchanged: 

N/2 



B = £ y/(k-l)(n + q)\n)(n-l | 
n=\ 

N-l 

+ £ y/(k-l)q + n(k+l)-N\n)(n-l \ + ^{k- \)q \ 0)(N- 1 | , 

n=f+l 

A = E l" + T>H- (4-38) 

In matrix representation: 

o o 




















































1 








R I \/(! + ?)(*- 1) " " i -j _ I u " " i (A K)l 

L > - ! y/(2 + q)(k-l) I' H ~ linn n l^- 1 ^ 



The corresponding density is shown in Fig. 7(a). 

A quasi-hole only affecting the upper layer of the v* = 2 fluid is shown in Fig.7(b). It is 
given by the solution: 

B = £ y/(k+l)(n+l) \n+l)(n\ + £ ^{k- \){a-q- 1 +nj \ n + l)(n | 

n=0 n=q+l 



iV-2 ^ 



£ ^(*+l)(p + n ___ (? _l)| n +l)< /I 



A = £l«)(?+l+«l+ £ | y+4+l+n)(24 + 2 + n| , (4.40) 
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Figure 7: Plot of the density of the v* = 2 Jain ground state, 1 /v d = k+ 1/2, for k = 4 and 
N = 400, including: (a) one quasi-hole in the origin (4.38) with q = 60; (b) a quasi-hole in 
the upper layer of the fluid (4.40) with q = 60; (c) the quasi-hole out of the origin (4.42) 
with q = 30 and r = 60. 
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with a = ^(z+g) ; p = 2+g ^ 1 ^ T and g a positive integer. In matrix representation for 
N = 8 and g = 1, it reads: 



/ 
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o 













y2(Fn) 
















/ 




1 \ 
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\ 1 / 
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/ 











(4.41) 



The displacement from the origin of the upper layer corresponds to Aj = (k + 1 ) (q + 1 ) . 

It is also possible to create a circular depletion in the whole fluid, of size (charge) Aj7 = q 
outside the origin at a distance Aj? = (k— l)(r— 1) (see fig. 7(c)): 



B 



A = 



r-l 



N ? 



£ y/(k-l){n + l) \n+l)(n\ + £ y/(k- l)(n+ 1 +q) \n + \){n 

n=0 n=r 



+ ^lj-q\r)(n-l\+\l^ + q)(k-l)\^)(^-l 



N-2 



+ £ ^(k+\)(n+\)-N+(k- \)q\ n+l)(n |, 



— ~i 



N . 

T 



£ \n+-)(n\. 

n=0 



(4.42) 



In matrix representation for N = 6 and q = 2 



B = 



A = 



( 



o 



o 
o 
o 
o 



o 
o 

o 










V(3 + r)(*-l) 





/ \ 





1 
1 

\ 1 / 








^J(A + r)k-2-r 




y/(5 + r)k-l-r 












(4.43) 



33 



In this case, the solution of the Gauss law is obtained in terms of a two-component auxiliary 
field \\f, and it holds for rq <C N, namely the depletion would move back to the origin in the 
scaling limit N — > °°. 

5 Conclusion 

In this paper, we have continued the study of the Maxwell-Chern-Simons matrix gauge 
theory as an effective theory of quantum Hall states. After providing better forms for the 
projection, A m 0, limiting state degeneracy, we have obtained the semiclassical ground 
states of the theory. They correspond to the quantum states found before [19], that repro- 
duce the Jain composite-fermion construction of phenomenological wave functions. The 
density of states in the main Jain series, V = m/ (mk+ 1), has been found to be that of in- 
compressible fluids: this confirms our expectation that the matrix states at g = are not too 
different from the physical states at g = °°. The semiclassical approximation used here is 
known to be valid in the large N limit for both matrix [22] [21] and real Hall states [29] [30] , 
in particular for incompressible fluid states. 

Outside the main series of hierarchical states, other ground states are possible in the ma- 
trix theory [19], that correspond to generalized Jain constructions [18]. In Jain's theory, 
these generalized states are excluded due to their low (or vanishing) gap. In the semiclas- 
sical analysis of the matrix theory, we have found that the majority of generalized states do 
not have piecewise constant density, i.e. are not incompressible fluids: this is an indication 
that they may become unstable for finite g > values. 

The study of the phase diagram of the matrix theory is clearly necessary to make better 
contact between the nice results (g = 0) and the physical regime (g = oo) 5 upon varying the 
potential V = g Tr[X,X] 2 . We expect that the relevant incompressible fluid states have a 
smooth evolution for g > and we plan to include the quartic potential in the semiclassical 
analysis by means of a mean-field approximation. 

The explicit semiclassical solutions in this paper can also be useful to study the symme- 
tries and algebraic properties of matrix ground states. We would like: 

• To make contact with the SU (m) symmetry of the conformal field theories describing 
the edge excitations of Jain states [31]. 

• To find a projection of states more refined than A m 0, that could discriminate the 
hierarchical Jain states from the generalized (unstable) ones. Such an expectation is 
based on the general belief that the observed Hall states should be uniquely character- 



34 



ized by algebraic conditions and gauge invariance, rather than by detailed dynamics, 
because they are exceptionally robust and universal. 
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A Appendix 

A.l Gauge invariance of the projection 

Here is an explicit proof that the projection (A a b) mx H (cf 2.24) is a gauge invariant condition 
on quantum states. Consider the more general relation for m = 2: 

A ab A^y (A,5) = M b y (A,fi) V a Wy . (A.l) 

The wave function is assumed to be gauge invariant: ^{A,B) = ^(UAU^ \UBU^). The 
form in the r.h.s. is specific of the bush states of section 2, but this is not relevant for 
the argument. The matrix Myy vanishes for a = a' \b = b' because *P is assumed to be 
one solution of the constraint. In general, there are several terms in the r.h.s. with that 
structure, but the matrices Myy should all vanish independently because they are multiplied 
by monomials V a Wy that are all independent [19]. 

Let us now multiply by unitary matrices and sum over indices on both sides to realize a 
gauge transformation of the two A's: 

{ UAU ') ab { UAU ") a >^^ = M rA^)ulut bi (uvuuw) a> , 

= M b y(lJAU\UBU^ (UV) a (UW) a , . (A.2) 

The resulting expression of M{UAJj\UBU^)by vanishes whenever M(A,B)t,y does, i.e. 
for b = b', because both vanish by polynomial identities that do not depend on the specific 
vales of the variables. Therefore, a solution of the constraint remains a solution after gauge 
transformation. 
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